Vortex Dynamics in a Coarsening Two Dimensional XY Model 
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The vortex velocity distribution function for a 2-dimensional coarsening non-conserved 0(2) 
time-dependent Ginzburg-Landau model is determined numerically and compared to theoretical 
predictions. In agreement with these predictions the distribution function scales with the average 
vortex speed which is inversely proportional to t^, where t is the time after the quench and x is near 
to 1/2. We find the entire curve, including a large speed algebraic tail, in good agreement with the 
theory. 



I. INTRODUCTION 



It is important to understand the role of defects in 
phase ordering [1] problems. We investigate here the 
growth kinetics of the non-conserved 0(2) symmetric 
time-dependent Ginzburg-Landau (TDGL) model in two 
dimensions after a quench from a disordered high tem- 
perature state to zero temperature. The dominant struc- 
tures in the ordering kinetics of this system are vortices 
with charges ±1. Vortices with higher order charges are 
unstable. Various aspects of the defect structure have 
been explored in some detail before [2,3]. We focus here 
on a numerical determination of the velocity distribution 
of the vortices as a function of time t after the quench. 
Theory predicts [4] that the distribution function scales 
with the average vortex speed which is inversely propor- 
tional to a length scale L{t), which grows with time t af- 
ter the quenches. One also finds a large speed algebraic 
tail in good agreement with predictions of an exponent 
of —3. In terms of the velocity distribution this corre- 
sponds to an exponent of —4. The number of vortices 
is also counted and its evolution in time is found to be 
consistent with previous work [5]. 

The disordering agents in the phase ordering of the 
n = d = 2 non-conserved TDGL model (2c? XY model) 
are well known, where d is the spatial dimension and n 
is the number of the components of the order parame- 
ter. One has unit one charged vortices and, at nonzero 
temperature, spin waves. If we focus on quenches to zero 
temperature we have only the vortices to consider. Thus 
a typical vortex configuration is shown in Fig. 1 . As the 
time evolves one has vortex anti-vortex annihilation un- 
til finally there are no surviving vortices and the system 
is fully ordered (we only consider the zero temperature 
case with no thermal fluctuations, where the system does 
eventually order). For general n = d Bray and Rutcn- 
berg [6] have shown that the growth law for such systems 
is given by L{t) « i^/^. The exception is for rt = c? = 2 
where their method is mute. The growth law for this case 
was treated by Pargellis et al. [7], and checked numeri- 
cally by Yurke et al. [5] . There is a logarithmic correction 
to the scaling law: L{t) « (i/ log(t))i/2. 
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FIG. 1. A typical vortex configuration in a 256 x 256 sys- 
tem with lattice spacing Ar = 7r/4. The arrow on each each 
site represents the order parameter at that point. Not all the 
lattice sites are shown. The squares and triangles are in the 
core regions of 4-1 and —1 vortices respectively, where the 
magnitude of the order parameter is near zero. The vortex 
core regions are picked out by using the method described in 
the text. 

There has also been theoretical work on the dynamics 
of these vortices. Mazenko [4] showed that if the order 
parameter can be assumed to be a gaussian field [8] when 
constrained to be near a vortex core then the vortex ve- 
locity probability distribution, for n = d [4,9], has the 
simple form: 



P{v) = 



r(l + n/2) 



(7rs2)«/2 (^_^^2/52)(«+2)/2 



(1) 



where the scaling speed s varies as for long times. 
If we only care about the magnitude of the velocity and 
integrate out the directions, then in the case of n = d = 2 
we have the speed probability distribution 



P(v) dv 
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or equivalently 



P{v) dv ■■ 
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(2) 



with a = (7r/2)^ and v = v/v, where the average speed 
V = TTs/2. So, after being scaled with the average speed, 
the vortex speeds have the same probability distribution 
at different times. A key feature of the predictions for 
P{v) is that there is a large velocity algebraic tail ^ v~^. 
This tail was also found using scaling arguments by Bray 
[10]. In this paper we check these predictions numerically 
for the case n = d = 2. We find that velocity probabil- 
ity distribution function does obey scaling of the form 
predicted by Eq. (2). But the average speed falls off as 
t~^^^ without the logarithmic correction and there is a 
large speed tail consistent with an exponent of —3. 

We have also monitored the vortex density Tly clS cl 
function of time and find, in agreement with previous 
work [5], 



(3) 



where Ly{t) oc (t/ log(t))^/^. The vortex number density 
should be proportional to the system's energy above its 
ground energy Eq (= — per unit area, where e is 
the control parameter of the system (see Eq. (4) below). 
Our simulation verifies this result. 



II. SYSTEM DESCRIPTION 

The system we study is described by the Langevin 
equation 



dipi 

'dt 



(4) 



where i = 1,2 are the indices for the two components of 
the order parameter ip. The noise term is zero because 
the system is quenched to zero temperature. By choos- 
ing proper units for the time and space and rescaling the 
order parameter, e can take on any positive value. The 
equation is put on a square lattice with periodic bound- 
ary conditions and driven by the finite difference scheme, 
i.e. replacing a* Vi(r,i) by [ip'fl^+^kl) - i^rikl)] /At with 
m being the time step number, and \/'^ipi{r,t) by 
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where A'^iV and NNN mean the nearest neighbors and 
next-nearest neighbors respectively and the lattice point 
is r/Ar = {k, I). Here At is the time step and Ar is the 
lattice spacing. Both are dimensionless. 

We have studied two systems in some detail. In both 
we choose e = 0.1 and use 1024 x 1024 lattice sites. In 



system one, or the bigger system, we use At = 0.02 and 
Ar = 7r/4. In system two, or the smaller system, we use 
At = 0.01 and Ar = tt/8. In both we measured the vor- 
tices number and the system energy. We measured the 
vortex speed distribution only in the bigger system. 

We prepare the system initially in a completely dis- 
ordered state. The average magnitude M of the vector 
order parameter ip at time t is calculated and the vortices' 
core regions are identified with those sites on which the 
order parameter magnitudes < M/4. Usually each 
core has about 10 sites in it. Here the coefhcient 1/4 is 
appropriately chosen so that no vortices are missed and 
no irrelevant points are picked. A circular integration 
around each vortex produces 27r or — 27r, which corre- 
sponds to two types of vortices of topological charges -1-1 
and —1 respectively. In fact there arc situations where 
we obtain charge 0. This is due to the non-zero area of 
the integration circle. When a pair of +1 and —1 vor- 
tices annihilate and the distance between them becomes 
smaller than the size of the integration circle, the circular 
integration will reflect the sum of the two charges, which 
is 0. Howevc^r, the lifetime of these charges are much 
smaller than the lifetime of the ±1 vortices. So they do 
not affect our statistics. We find in our simulations that 
the numbers of positive and negative vortices are equal. 

The position of a vortex is given by the center of its 
core region. Suppose the order parameter's magnitude 

at the core region is described by M{xi,yi) with {xi,yi) 
belonging to the core region. Then by fitting M{xi,yi) 
to the function M{x, y) = A + B[{x - xof' + [y - yo)^] 
we can find the center (xo,yo)- 

The positions of each vortex at different times are 
recorded, and the speed is calculated using v = Ad/ At. 
Here Ad is the distance that the vortex travels in time 
Ar. We have, simultaneously, recorded the number of 
vortices as a function of time and checked the scaling 
result given by Eq. (3). 



III. NUMERICAL RESULTS 



Every 5t = 10 we compute the speed of each vortex 
with At = 5, i.e. v = \r{t + At) - r(i)|/AT with t in- 
creasing by step length 5t and r being the position of 
the vortex. We found the average speed of the vortices 
v{t) is proportional to t~^ with x = 0.51 ± 0.01 for the 
bigger system (Ar = 7r/4) as shown in Fig. 2. Unlike 
the case for Ly{t) extracted from Uy and {E — Eq), we 
do not observe logarithmic correction for v{t). We can 
not rule out such corrections appearing on a longer time 
scale but they do not enter on the same time scale as for 

Ly{t). 
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FIG. 2. The average speed v{t) of the vortices for the bigger 
system (Ar — 7r/4). The speed of each vortex is calculated 
with Ar = 5. t"^ is used to fit the data, x = 0.51 ± 0.01 for 
the system with Ar = n/4. The data are averaged over 60 
different initial conditions. 
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FIG. 3. The vortices' speed distribution probability den- 
sity for the Ar = 7r/4 system P{v/v(t)) (the solid line) fit 
to the function Pix) — 2ax/{l + aa;^)^ with a = 2.12 and 
X = v/v{t). The theoretical curve with a — (it/ 2)^ is also 
shown in the figure. The large speed tail of the distribution 
can be fit to with y = 3. The insert shows the same data 
in the logarithmic scale, and it seems that the the theoretical 
curve fits the tail better. The data are averaged over 60 dif- 
ferent initial conditions. In the calculation of the probability 
distribution we use the bin width 0.02. 



We compared the speed distribution at different times. 
They have approximately the same shape after being 
rescaled by the average speed. At early times this is 
clear. At late times there are fewer data points and the 
similarity between the two distributions at different times 
is not so apparent. Even at early times the number of 
data points are not enough to give a good fit to the dis- 
tribution's long tail. So we rescale the speed data with 
the best fit to the average speed v{t) oc t~^-^^ and put the 
data for all times into one histogram as shown in Fig. 3. 
By this means we obtain better statistics. We find that 
the distribution can be well fit to the function 



In these results we have scaled the velocity with ^ 
by taking z = x just as the theory predicts. However 
the value of scaling exponent z of the speed distribution 
is not very robust in our simulations. If we change the 
exponent z and use t^^ to rescale the speed distribution 
at different times, and they also have approximately the 
same shape. So the uncertainty of z is quite large. After 
being rescaled by t^^ and put into one histogram, the 
speed distribution for small and intermediate values of 
V are somewhat dependent on the value of z. However, 
the power-law tail is quite robust. We alter the value of 
z between 0.4 and 0.6, and find that the tail exponent 
changes by less than 0.1. 



P{v) 
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(6) 



for a between 2 and 2.5. The best fit is for a = 2.12. How- 
ever if we require that the average speed is given by the 
measured speed then we must have a = (7r/2)^ = 2.47. 
The best fit and the most consistent fit are both shown 
in Fig. 3. The distribution has a long tail which is ap- 
proximately [v /v)~'^ with y = 3. If wc just try to fit the 
long tail, then the best fit is given by a = 2.47. These 
results are in excellent quantitative agreement with the 
theoretical prediction in Ref [4] . 
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FIG. 4. The total number of the vortices per unit area 
nv(t)/S times t, where S — (1024 x Ar)^. In both systems 
Hy ■ t/S can be fit to alog{t/to)- The data from the laxge and 
the small systems are averaged over 68 and 58 different initial 
conditions respectively. 
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FIG. 5. The product of the systems energy per unit area 
{E — Eo)/S with t. The ground state energy is Eq = —e^S/A. 
The data can also be fit to a ln(t/io). The data from the large 
and the small systems are averaged over 68 and 58 different 
initial conditions respectively. 

In both the larger and smaller systems the vortex num- 
ber densities have the same time dependence. We obtain 
riy ~ [i/ log(t/to)]~"^, where is the number density of 
vortices (positive or negative) and to is a constant. In 
Fig. 4, we show the data for n,, • t/S. S is the area of 
the system. In Fig. 5, we show the data for the energy 
per unit area {E — Eq)/S times t. We conclude that the 
energy density is proportional to the number density of 
the vortices. 

IV. CONCLUSIONS 

We have studied the growth kinetics of the non- 
conserved TDGL model in the case of n = d = 2. We 
measured the speed probability distribution for the ±1 
vortices. At any given time with relatively few vortices, 
the statistics are poor. However the accumulated data for 
all times when scaled gives a scaling function with good 
statistics. Although the scaling exponent z is 0.5 with 
significant uncertainty, the large speed tail does takes 



the form of (v/v)'^ with the exponent y = 3. This is 
consistent with the theoretical prediction. The form of 
the distribution P{v) is quantitatively consistent with the 
theoretical prediction. 

Why does the theory do so well? It was shown by 
Mazenko and Wickham [11], that one can construct a 
nontrivial self- consistent gaussian theory for the order 
parameter if it is constrained to be evaluated near a vor- 
tex core. Such constraints occur naturally, for example, 
in averages over the vortex density. This suggests that 
the theory developed in Ref. [4] may be on a firmer foot- 
ing than first thought. 

According to the theory, the average speed should be 
proportional to the inverse of the correlation length L{t). 
We did not observe any logarithmic correction in the av- 
erage speed, although it appears in the correlation length 
when n = d = 2. However a logarithmic behavior may 
still exist. The time we used in our simulations may not 
be long enough to see the effect. 

We also measured the number density of the vortices 
and the energy density of the two systems. Their time 
dependences both have a logarithmic correction, which is 
consistent with previous work. 
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